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ABSTRACT 

In the present work, we investigate the universe dominated by quintom or hessence energies in 
Loop Quantum Cosmology (LQC). Interestingly enough, we find that there are some stable attrac- 
tors in these two cases. In the case of quintom, all stable attractors have the feature of decelerated 
expansion. In the case of hessence, most of stable attractors have the feature of decelerated expan- 
sion while one stable attractor can have decelerated or accelerated expansion depend on the model 
parameter. In all cases, the equation-of-state parameter (EoS) of all stable attractors are larger than 
— 1 and there is no singularity in the finite future. These results are different from the dynamics of 
fT^ ■ phantom in LQC, or the ones of phantom, quintom and hessence in classical Einstein gravity. 
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I. INTRODUCTION 

Dark energy fl*] has been one of the most active fields in modern cosmology since the discovery of 
accelerated expansion of our universe 0, 0, d, 0; 13 ■ the observational cosmology of dark energy, 
equation-of-state parameter (EoS) Wde = Pde/pde plays an important role, where p^e and pde are the 
pressure and energy density of dark energy respectively. Recently, evidence for Wde iz) < —1 at redshift 
z < 0.2 ~ 0.3 has been found by fitting observational data (see [1 [13, [H [ll [ll, Q d, [ll] for 
examples). In addition, many best-fits of the present value of Wde are less than —1 in various data fittings 
with different parameterizations (see for a recent review). The present data seem to slightly favor an 
evolving dark energy with Wde being below —1 around present epoch from Wde > — 1 in the near past [13] ■ 
Obviously, the EoS cannot cross the so-called phantom divide Wde = — 1 for quintessence or phantom 
alone. Some efforts have been made to build dark energy model whose EoS can cross t he p hantom divide 
(see for examples [13, [ll [IS [13, [11 [11, [13, [13, [11, [IS [13, [11, [H, M, M, M, M, M, [H] and references 
therein). 

In [lO|, Feng, Wang and Zhang proposed a so-called quintom model which is a hybrid of quintessence 
and phantom (thus the name quintom). It is one of the simplest modeles whose EoS can cross the phantom 
divide. The cosmological evolution of the quintom dark energy was studied in [13, [13| . Perturbations 
of the quintom dark energy were investigated in [^56, 21}', and it is found that the quintom model is 
stable when EoS crosses —1, in contrast to many dark energy models whose EoS can cross the phantom 
divide [11]. Other works concerning quintom also include |3l| for examples. In [13], by a new view of 
quintom dark energy, one of us (H.W.) and his collaborators proposed a novel non-canonical complex 
scalar field, which was named "hessence", to play the role of quintom. In the hessence model, the 
phantom-like role is played by the so-called internal motion 0, where 9 is the internal degree of freedom of 
hessence. The transition from ii;>— ltow<— lor vice versa is also possible in the hessence model (isj . 
The cosmological evolution of the hessence dark energy was studied in [Tol] and then was extended to the 
more general cases in [23|. The w-w' analysis of hessence dark energy was performed in [2l| . In [33], 
the method to reconstruct hessence dark energy was proposed. We will briefly review the main points of 
quintom and hessence energies in Sec. [iTl 

In fact, many works by now are considered in the framework of classical Einstein gravity. However, it is 
commonly believed that gravity should also be quantized, like other fundamental forces. As well-known, 
for many years, the string theory is the only promising candidate for quantum gravity. In the recent 
decade, however, the Loop Quantum Gravity (LQG) (see e.g. [H, [13, [13, [lH for reviews) has became a 
competitive alternative to the string theory. LQG is a leading background independent, non-perturbative 
approach to quantum gravity. At the quantum level, the classical spacetime continuum is replaced by 
a discrete quantum geometry and the operators corresponding to geometrical quantities have discrete 
eigenvalues. 

Loop Quantum Cosmology (LQG) (see e.g. [13, [13, [11] for reviews) restricts the analysis of LQG to 
the homogeneous and isotropic spacetimes. Recent investigations have shown that the discrete quantum 
dynamics can be very well approximated by an effective modified Friedmann dynamics [45', '46', [471 . 
There are two types of modification to the Friedmann equation due to loop quantum effects [44, ]54| . 
The first one is based on the modification to the behavior of inverse scale factor operator below a 
critical scale factor a*. So far, most of the LQG literature has used this one. Many interesting results 
have been found, for instance, the replacement of the classical big bang by a quantum bounce with 
desirable features [13, [13, [13] , avoidance of many singularities [13], easier inflation [sol. [5l| . correspondence 
between LQG and braneworld cosmology [HI], and so on. However, as shown in e.g. [53j]. the first type 
of modification to Friedmann equation suffers from gauge dependence which can not be cured and thus 
lead to unphysical effects. In a recent paper [63], Magueijo and Singh provided very sharp results to show 
that it is only for the closed model that such modifications can be sensible and for flat models they make 
no sense. 

The second type of modification to Friedmann equation is discovered very recently [i^ . [ssl . [s^ , which 
essentially encodes the discrete quantum geometric nature of spacetime. The corresponding effective 
modified Friedmann equation in a flat universe is given by [H, [53, [sJ, [s^, [sl, [53, [13] 
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where H = a/ a is the Hubble parameter; a dot denotes the derivative with respect to cosmic time i, and 
a is the scale factor; p is the total energy density; = SttG; and the critical density reads 

^ (2) 



167rVG2;i 



where 7 is the dimensionless Barbero-Immirzi parameter (it is suggested that 7 ~ 0.2375 by the black 
hole thermodynamics in LQG ^56*1). Differentiating Eq. H]) and using the energy conservation equation 

p + iH[p + p) = 0, (3) 

one obtain the effective modified Raychaudhuri equation [s^, [s^, [H, [H^ [lOl 

ij^-^(p + p)fl-2^V (4) 



2 " ^ ' V A 

where p is the total pressure. Actually, as shown in [63 | , the effective modified Raychaudhuri equation can 
be derived directly by using the Hamilton's equations in LQC, without assuming the energy conservation 
equation. It is easy to check that only two of Eqs. H]), ^ and @ are independent of each other, and 
the third one can be derived from these two. By using the second type of modification to Friedmann 
equation, the physically appealing features of the first type are retained 14^ 1 for instance, resolution of 
big bang singularity avoidance of big rip and other singularities [S^jSjIs^j inflation in LQC |,60| . 

correspondence between LQC and braneworld cosmology [54[ , and so on. 

For the universe with a large scale factor, the first type of modification to the effective Friedmann 
equation can be neglected and only the second type of modification is important (53 | , while the matter 
Hamiltonian T-Lm and the corres pon ding expressions for energy density and pressure retain the same 
classical forms [i^. Issl. Is^. Issl IstI. ISSl [59l| . This is the particular case which we will consider here. 

In the present work, we will investigate the universe dominated by quintom or hessence in LQC. 
Following [53, [13, m, , we use the method of dynamical system foTl . After a brief review of quintom 
and hessence energies in SecHH we consider the dynamics of quintom and hessence in Sec. lIIII and Sec.HVl 
respectively. Interestingly enough, we find that there are some stable attractors in these two cases. In the 
case of quintom, all stable attractors have the feature of decelerated expansion. In the case of hessence, 
most of stable attractors have the feature of decelerated expansion while one stable attractor can have 
decelerated or accelerated expansion depend on the model parameter. In all cases, the EoS of all stable 
attractors are larger than —1 and there is no singularity in the finite future. These results are different 
from the dynamics o f ph antom in LQC (58l.[59l|. or the ones of phantom, quintom and hessence in classical 
Einstein gravity [ll, [20|, jH, [H HI] . 



II. QUINTOM AND HESSENCE ENERGIES 



A. Quintom energy 

Phenomenologically, one may consider the Lagrangian density for quintom [lol . [23l [23 | 

A = ^ (<9^'/'l)' - \ (5p02)' - ^(01>2), (5) 

where ^1 and 02 are two real scalar fields and play the roles of quintessence and phantom respectively. 
Considering a spatially flat Friedmann-Robertson- Walker (FRW) universe and assuming the scalar fields 
01 and 02 are homogeneous, one obtains the effective pressure and energy density for the quintom, i.e. 

P=\4>\-\^l-V{(j^lA2). p=i0?-i02 + y(0i,02), (6) 
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respectively. The corresponding effective EoS is given by 

P 4>l-4>l + 2V{(t)i,(l)2)' 



(7) 



It is easy to see that w > —1 when > (/)| while w < —1 when < cp^- We consider the simplest 
quintom whose V{4>1t4'2) = Vi{cj)i) + V2(02) in the present work. In this case, the equations of motion 
for 01 and 02 are given by 

01 + 3i/0i + ^ = 0, (8) 

02 + 3i/02 - ^ = 0. (9) 

d02 



B. Hessence energy 

Following [H, [l^ [lOl , we consider a non-canonical complex scalar field, namely the hessence, 

$ = 01+i02, (10) 

with a Lagrangian density 

Ch^\[ id,<i>r + id.^r ] ui^' + 'f*') = \ [ id.cbf - cj^^d.er ] - ^(0), (n) 

where we have introduced two new variables (0, 9) to describe the hessence, i.e. 

01 = 0cosh0, 02 = 0sinh^, (12) 

which are defined by 

bl-(j)l, coth6i = ^. (13) 



a2 



In fact, it is easy to see that in terms of 0i and 02, the hessence can be regarded as a special case of 
quintom with general y(0i,02). Considering a spatially flat FRW universe with scale factor a{t) and 
assuming and 9 are homogeneous, from Eq. (fTTj) we obtain the equations of motion for and 9, 

+ 3i70 + 002 + ^ = 0, (14) 

rf0 

0^0 + (200 + 3i/02)^ = 0. (15) 
The pressure and energy density of the hessence are 



2 V" " / 2 
respectively. Eq. ([T5| implies 

Q = a^ct)'^9 = const. (17) 

which is associated with the total conserved charge within the physical volume due to the internal 
symmetry [l^, [l^ . It turns out that 



a-^0^ 
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Substituting into Eqs. p4|) and ^6]). they can be rewritten as 

Noting that the EoS w = p/p, it is easy to see that w > —1 when (fP > Q"^ / [a^ cfp') , while w < —1 when 
</)^ < Q'^ /{a^cj}^). We refer to the original papers [lljEllSl for more details. 



III. DYNAMICS OF QUINTOM ENERGY IN LQC 

In this section, we consider the universe dominated by quintom energy in LQC. Following (23l. [23. Issl. 
[59L [5^ , we introduce these five dimensionless variables 

K(j)l K<j)2 K^/V\ KyfSAl p 

We introduce z just for convenience. It is expected that z is not independent, because of Eq. In 
fact, in the case of the universe dominated by quintom energy in LQC, the effective modified Friedmann 
equation, namely Eq. ([1]), can be rewritten as 

{xl-xl+yl + yl){l-z) = l. (22) 

Thus, one can explicitly express z in terms of xi, X2, yi and y2- Notice that < z < 1 is required by the 
positiveness of p and in Eq. ([T]). In addition, by using Eq. ([S]), we recast Eq. ^ as 



= 3 (x? - xl) (1 - 2z) , (23) 



which will be used extensively. Also, Eq. ([7]) becomes 

xl-xl+yl+vl 

In this work, we consider the case of quintom with exponential potentials 

VM^) = V^,e-^^''^\ ^2(02) = V^0.e-^^"*% (25) 

where Ai and A2 are dimensionless constants. Without loss of generality, we choose Ai and A2 to be 
positive, since we can make them positive through field redefinition 0i — > —(pi^ (t>2 —4>2 if Ai and A2 
are negative. By the help of Eqs. ([22|) — ((25|) . ^ and ([3]), the equations of motion for 0i and (t>2, namely 
Eqs. ([5]) and ([5]), can be rewritten as an autonomous system 

x'l = 3x1 (^ - 1) + (26) 

x'2^ix2{A-l)- ^j^\2yl, (27) 

y[ = iy,I^A~^x,y (28) 

J/2 = 32/2 [a - ^x^ , (29) 

z' = ~Qz{x\-xl){l- z), (30) 
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where a prime denotes the derivative with respect to the so-caUed e-folding time iV = In a, and 



A: 



in which we have used Eq. 



\){\~2z) = {xl~xl) 2{xl 



vl 



2/2 



(31) 



We can obtain the critical points {xi,X2,yi,y2T z) of the autonomous 
system Eqs. p6|) — (pO| by imposing the conditions x'l = x'2 = y[ = y'2 = z' = 0. Of course, they are 
subject to the Friedmann constraint Eq. (f^ . namely (a;^ — X2 + yf + yf) (1 — ^) = 1- We present the 
critical points and their existence conditions in Table [H In fact, there are other four critical points with 



Xi 



A1A2 



X2 



XIX2 



yi ^ ±A2 



'9 ' 



y2 = ±Xl^/^1 



q 1 



z = 0, 



where 



%/6(Ai-A?)' 

^ 6Ai-A^ (6 + Ai) 

Hxi-xlf 

However, the existence of yi and y2 requires = 0. In this case, these four critical points reduce to 
Points (Q2p) or (Q4p). 



Label 


Critical Point {x\,X2,y\,y2, z) 


Existence 


Qip 


x'i > 1, ^xf - 1, 0, 0, 


xl>l 


Qlm 


Si > 1, -\/s? - 1, 0, 0, 


xl>l 


Q2p 
Q2m 


f , y^l-i, 0, 0, 
-v^i|-i. 0, 0, 


Xi<VE 


Q3 
Q4p 


^. 0, ±\/l-#, 0, 
71 + 4, ^, 0, 0, 


Xi<V6 
always 


Q4m 




always 


Q5 




always 



TABLE I: Critical points for the autonomous system Eqs. I|26|) — (|30p and their existence conditions 



To study the stability of the critical points for the autonomous system Eqs. ()26p — (PU)) . we substitute 
linear perturbations — > xi + 6x1, X2 ^ X2 + Sx2, yi ^ yi + Syi, y2 ^ y2 + 5y2, and z ^ z + 5z 
about the critical point (xi, a;2, yi, ^2, z) into the autonomous system Eqs. (|26p — pop and linearize them. 
Because of the Friedmann constraint (I22p . there are only four independent evolution equations, i.e. 

(32) 
(33) 

(34) 



5x'^ = ixiSA + 3 (A - 1) 5xi + V&Xiyi5yi, 
5x2 = 3X2M + 3 (A - 1) 5x2 - V6A2272<52/2, 



where 



6yi = 3?;i ( 5 A - ^^^1 



5y2 = 3?;2 5 A 



5x' 



Syi, 
Sy2, 



(35) 



A -~ 
5A 



yxi X2) 
-4 {xl - a 
2{xl 



2{xl 



yl 



vl) 



vl) 



vl) 



- 1 



{xi5xi — X26X2 + yi5yi 

{xiSxi - X25X2) ■ 



V25y2) 



(36) 



(37) 
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The four eigenvalues of the coefhcient matrix of Eqs. (|32|) — ([35|) determine the stabihty of the critical 
point. We present the corresponding eigenvalues for the critical points in Table [TTl It is easy to see 
that Points (Qlm), (Q2m), (Q3), (Q4m) and (Q5) are unstable. Point (Qlp) exists and is stable under 
conditions xi > 1, Ai^i > and A2\/sf — 1 > a/G. Point (Q2p) exists and is stable under conditions 



Ai < V6 and X2y^^i ^ — 3/2 > 3. Point (Q4p) exists and is stable under condition Aiy SAj ^ + 3/2 > 3. 

The three stable attractors (Qlp), (Q2p) and (Q4p) have the common features yi — y2 = z — and 
x\ — x^ — 1. From Eq. ([24|) . we get the EoS w = 1, which implies that the quintom behave as a stiff fluid. 
From Eq. ([^5]) . ~H/H^ — 3. Then, we find that H = /3 (the integral constant can be set to zero by 



redefining the time). Thus, we obtain a oc t^^^. From Eq. ^ and w — 1, we find that p 
The universe undergoes decelerated expansion and there is no singularity in the finite future. 



oc t 



Point 



Eigenvalues 



3- y^fAiii, 3-A2^/|(a;2-l) 



Qlp 

Qlm 

Q2p 

Q2m 

Q3 

Q4p 

Q4m 

Q5 



-6, 
-6, 0, 3 



3 + Aay^ 



-6, 0, 0, 3 -A; 
-6, 0, 0, 3 + A2y^ 
-A?, I (Af - 6), i (A 



-6, 0, 0, 3-Aiy'9A 
-6, 0, 0, 3 + Aiy9A 



-2,3 

2 2 

-2,3 

2 2 



Ai, -H^2 + 6), -i(Ai+6) 



TABLE II: The corresponding eigenvalues for the critical points of the autonomous system Eqs. (|26p — (|30p . 



IV. DYNAMICS OF HESSENCE ENERGY IN LQC 

In this section, we consider the universe dominated by hessence energy in LQC. Similar to the case of 
quintom, following [H, [13, [H, [sl, [111 , we introduce these five dimensionless variables 

k4> kVV p \/6 K, Q 

y=^?77' z=—, u=—, z;=^=— — . (38) 



Again, we introduce z just for convenience, since it is expected that z is not independent due to Eq. (|20p . 
In fact, in the case of the universe dominated by hessence energy in LQC, the effective modified Friedmann 
equation, namely Eq. ([T]), can be rewritten as 

{x'' -v''){l-z)^l, (39) 

which can be used to explicitly express z in terms of x, y and v. Notice that < z < 1 is required by the 
positiveness of p and H"^ in Eq. By using Eq. (^0]) . the effective modified Raychaudhuri equation, 
namely Eq. can be rewritten as 

-^ = i{x^^v^){l-2z). (40) 
From Eq. ([20)1 . the EoS of hessence is given by 

7 9 7 

p x^-v^-y^ 
p x^ ~ + y^ 
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In this work, we consider the case of hessence with exponential potential 



(42) 



where A is a dimensionless constant. Without loss of generality, we choose A to be positive, since we can 
make it positive through field redefinition (f) —<j) if A is negative. By the help of Eqs. ((39|) — (I42p . (f20|) 
and ([3]), the equation of motion for 0, namely Eq. can be rewritten as an autonomous system 



y' ~ 3y ( B —x 

z' — — 6z (a;^ — w^) 
u' = —xu^, 



il-z) 



3v \ B — 1 — 



where 



B={x^- v^) (1 - 2z) = (a 



2 x' 



(43) 

(44) 

(45) 
(46) 

(47) 
(48) 



in which we have used Eq. (j39p . We can obtain the critical points (x, y, z, u, v) of the autonomous system 



Eqs. — (|T7|) by imposing the conditions x' 
the Friedmann constraint Eq. ([39|) . nan 
their existence conditions in Table [TTTl 



= z' = u' = v' = 0. Of course, they are subject to 
w^) {I — z) = 1. We present the critical points and 



Label 



HI 

H2p 

H2m 

H3 

H4 



Critical Point {x, y, z, u, v) 



x'^ > 1, 0, 0, 0, ±V2;2 - 1 
1, 0, 0, 0, 



V6 



, 0, 0, 0, ±/ 



^, 0, 0, 



Existence 



> 1 
always 
always 

x<Vg 

A < \/6 



TABLE IIL Critical points for the autonomous system Eqs. (|43p — H47p and their existence conditions. 



To study the stability of the critical points for the autonomous system Eqs. (|43p — (|T7|. we substitute 
linear perturbations a; — > a; + (5a;, y y + Sy, z z + 6z, u — > u + Su, and d ^ ti + (5v about the 
critical point {x,y,z^u,v) into the autonomous system Eqs. (|43p — (|47p and linearize them. Because of 
the Friedmann constraint (j39p . there are only four independent evolution equations, i.e. 

(49) 
(50) 
(51) 
(52) 



Sx' — 3xSB + 3 (-B — l) (5x — 2uvSv — v^Su + VGXySy, 
Sy' - 3y (^6B - -^Sx^ + 3 (^B - -^x^ Sy, 



6u' 

Sv' — 3v 



2xu6u — u^Sx 
1 



SB — — {xSu + uSx) 



3 LB- 1 



1 



xu ] Sv, 



where 



5= (x^-v^) [2 {x^ + y^ - v^' - l\ , 

SB = —4 (x^ — v^) (x^ + j/^ ^ v'^) ^ {xSx + ySy — vSv) 
7^) "'^ — 1 {xSx — vSv) . 



^2 2(x^+f 



(53) 
(54) 
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The four eigenvalues of the coefhcient matrix of Eqs. (|49p — ([52|) determine the stabiHty of the critical 
point. We present the corresponding eigenvalues for the critical points in Table [iVl It is easy to see 
that Point (HI) exists and is stable under conditions x > 1 and Xx > Point (H2p) is stable under 
condition A > \/6; Point (H2m) is unstable; Points {H3) and (H4) are stable when they exist under 
condition A < \/6- 



Point 


Eigenvalues 


HI 


-6, 0, 0, 




H2p 


-6, 0, 0, 3 - J 




H2m 


-6, 0, 0, 3+^ 


|A 


H3 


-6, 0, 0, 




H4 


0, -A^ H^'-6), i 


(A^ - 6) 



TABLE IV: The corresponding eigenvalues for the critical points of the autonomous system Eqs. (|43p — (|47|l . 

The stable attractors (HI), (H2p) and (H3) have the common features y = z = u — and — = 1. 
From Eq. (|4T|) . it is easy to find that the EoS w = 1, which implies that the hessence behave as a stiff 
fluid. From Eq. ([3D]), —H/H^ = 3. Then, we find that H = /3 (the integral constant can be set to zero 
by redefining the time). Thus, we obtain a ex t^^^. From Eq. ^ and w = 1, we find that p cx cx t~^. 
The universe undergoes decelerated expansion and there is no singularity in the finite future. 

The stable attractor (H4) is slightly different from other three stable attractors. From Eq. ([4T|) . we 
get the EoS w = -1 + X^/S > -1. Note that w > -1/3 for V2 < A < V6, while w < -1/3 for A < ^2. 
From Eq. ~H/H^ = A^/2. Then, wc find that H = 2t^^ j)? (the integral constant can be set to 

zero by redefining the time). Thus, we obtain a oc t^^^ . From Eq. ([3]) and w — —1 + A^/3, we find 
that p cx oc t~^. The universe experiences decelerated expansion for < A < VG, or accelerated 
expansion for A < \/2. However, the universe cannot undergo super-accelerated expansion [H > 0) for 
any A. Therefore, there is no singularity in the finite future for any case. 



V. CONCLUSION 

In the framework of classical Einstein gravity, the dynamics of phantom, quintom and hessence have 
been studied in literature [H, HO, IH, IH, Ml] . In the case of phantom, the universe will end in a big rip 
singularity [6^ . In the case of quintom without direct couping between (f>i and (j)2 11^ , or with a special 
interaction between 0i and (j)2 through Vint [Vi {(jji) V2 {4>2)]^^'^ [111, the phantom-dominated solution 
is the unique attractor and the big rip is inevitable. In the case of hessence [H, [l^], however, the big rip 
can be avoided, while its attractors are slightly different from the ones of the present work. 

In the framework of LQC, the dynamics of phantom has also been studied [58, 59]. It is found that there 
is no stable attractor in this case. Therefore, the phase trajectory is very sensitive to initial conditions. 
However, the big rip can be avoided and the universe finally enters oscillatory regime. This is mainly due 
to the quantum correction to Friedmann equation. 

In the present work, we investigate the universe dominated by quintom or hessence energies in Loop 
Quantum Cosmology (LQC). Interestingly enough, we find that there are some stable attractors in these 
two cases. In the case of quintom, all stable attractors have the feature of decelerated expansion. In 
the case of hessence, most of stable attractors have the feature of decelerated expansion while one stable 
attractor can have decelerated or accelerated expansion depend on the model parameter. In all cases, the 
equation-of-state parameter (EoS) of all stable attractors are larger than —1 and there is no singularity 
in the finite future. These results are different from the dynamics of phantom in LQC, or the ones of 
phantom, quintom and hessence in classical Einstein gravity. 
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